function r = myh(h,k,x)
dxdx=2*k+8*h*x(1)^2+4*h*(x(1)^2+x(2)-11);
dxdy=4*h*x(1)+4*k*x(2);
dydx=dxdy;
dydy=2*h+8*k*x(2)^2+4*k*(x(1)+x(2)^2-7);
r = [dxdx, dxdy;
     dydx, dydy ];
end